Intramolecular vibrational energy redistribution in DCO (X 2 A'): Classical-Quantum 
correspondence, dynamical assignments of highly excited states, and phase space 

transport 
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Intramolecular dynamics of highly excited DCO (X 2 A') is studied from a classical-quantum cor- 
respondence perspective using the effective spectroscopic Hamiltonian proposed recently by Trollsch 
and Temps (Z. Phys. Chem. 215, 207 (2001)). This work focuses on the polyads P = 3 and P = 4 
corresponding to excitation energies E v ~ 5100 cm -1 and 7000 cm -1 respectively. The majority 
of states belonging to these polyads are dynamically assigned, despite extensive stochasticity in the 
classical phase space, using the recently proposed technique of level velocities. A wavelet based 
time- frequency analysis is used to reveal the nature of phase space transport and the relevant dy- 
namical bottlenecks. The local frequency analysis clearly illustrates the existence of mode-specific 
IVR dynamics i.e., differing nature of the IVR dynamics ensuing from CO stretch and the DCO 
bend bright states. In addition the role of the weak Fermi resonance involving the CO stretch and 
DCO bend modes is investigated. A key feature of the present work is that the techniques utilized 
for the analysis i.e., parametric variations and local frequency analysis are not limited by the di- 
mensionality of the system. This study, thus, explores the potential for understanding IVR in larger 
molecules from both time domain and frequency domain perspectives. 



I. INTRODUCTION 



Substantial progress has been made over the past few 
decades in understanding the nature of energy flow in 
highly excited molecule»i*2i 3 ^i£*iii&. At higher levels of 
excitation the usual rigid rotor-normal mode descrip- 
tion breaks down and the modes get strongly coupled. 
The resulting spectral complexity in terms of the split- 
tings and intensities encode the flow of energy through 
the molecule. This phenomenon of intramolecular vibra- 
tional energy redistribution (IVR) can lead to subtle but 
important deviations from the predictions of statistical 
theories like RRKMiS and thus profound consequences 
for the rates of chemical reactions 9 . The mode mix- 
ings in general have classical as well as quantum origins. 
Much of the classical effects arise from vibrational an- 
harmonicities, coriolis and centrifugal distortion effects. 
Possible quantum contributionsii*i2iiii^iiSii& manifest as 
indirect state-to-state explorations involving a sequence 
of intermediate, off-resonance virtual states ("vibrational 
superexchange" ) . In molecules with symmetry the phe- 
nomenon of dynamical tunneling provides a quantum 
route to mode mixingiii^i&. Theoretically it is cru- 
cial to establish and understand IVR pathways from 
time dependent as well as time independent perspectives. 
RRKM theory has been quite successful^ but evidence is 
mounting for a number of reactions where the intramolec- 
ular dynamics is inherently nonstatistical 9-10 . Account- 
ing for such nonstatistical dynamics is theoretically chal- 
lenging and crucial for control of reactions. 

In this work our aim is to understand IVR in the 
deuterated formyl radical DCO. Our choice of the system 
is primarily due to the considerable experimental^ and 
theoreticali£ii2i22i2i work on IVR and mode specificity 
of the unimolecular decay of highly vibrationally excited 



DCO. In contrast to HCO the deuterated analog involves 
a strong accidental 1:1:2 resonance between the vi- 
brational frequencies. The experiments and analysis of 
Renth et al using a effective Hamiltonian clearly shows 
that the unimolecular decay in DCO is IVR limited^. 
Although DCO has three vibrational modes, m (DC 
stretch), n (CO stretch), b (bend), the coupling structure 
of the effective Hamiltonian^ yields a conserved quantity 
called as the polyad number P = m + n + b/2. Ex- 
perimentally the existence of P is revealed by bunching 
of vibrational levels as evident in the SEP spectrum of 
DCOii. 

As a small molecule exhibiting nonstatistical in- 
tramolecular dynamics DCO provides an ideal system 
to test the various theoretical approaches for study- 
ing highly excited molecules. Impressive advances have 
been made towards fully quantum investigation of the 
IVR dynamics in large systems£i 9 ^i2i2^. There is little 
doubt that DCO can be subjected to such exact treat- 
ments which are important in their own right. How- 
ever there are strong reasons^ to believe that pow- 
erful insights into IVR can be gained by approaches 
based on classical-quantum correspondence. For exam- 
ple, Jung et al have provided dynamical assignments 
of the highly excited states of DCO based on certain 
key periodic orbits in phase spaceSI. A key to the 
success of the many classical-quantum correspondence 
studiea 26 . 27 i 28 i 29 i 30 . 31 i 32 i 33 i 34 i 35 i 36 is that the understand- 
ing of IVR is firmly rooted in the invariant structures 
inhabiting the underlying phase space. This provides a 
very clear dynamical picture of IVR in molecules and can 
serve as a baseline for estimating the important quantum 
routes to mode mixings. 

Currently one of the main obstacles in extending such 
useful classical-quantum correspondence studies is the is- 
sue of dimesionality. Despite the existence of potentially 
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3iV — 6 coupled vibrational modes it is now well under- 
stood that not all of the modes are coupled strongly^ 9 .. 
Although such hierarchicality of the couplings leads to 
much smaller dimensionalities of the 'IVR space' it is pos- 
sible that for large enough molecules or smaller molecules 
at higher levels of excitation the effective dimesional- 
ity of the system can be larger than three. In such 
instances the dimensionality of the system poses both 
technical as well as conceptual problems in the context 
of classical-quantum correspondence studies. The tech- 
nical problems are well known and have to do with the 
inability to visualize Poincare surface of sections or the 
quantum states and, computation of the periodic orbits 
or higher dimensional tori. Thus the knowledge of the 
global phase space structure along with its influence on 
the quantum states is at present difficult to obtain for sys- 
tems with three or more degrees of freedom. For more 
detailed discussions on the dimensionality issues, tech- 
nical and conceptual, we refer the reader to the earlier 
-.-□rfa 29 - 37 - 38 - 39 ' 40 . 

One of the central goals of our research is to come 
up with techniques that would allow a detailed classical- 
quantum correspondence study of IVR in molecules with 
three or greater degrees of freedom and in this context 
two approaches will be of particular interest. The first 
one is the promising technique of local frequency analysis 
(LFA) put forward nearly two decades ago by Martens, 
Davis and Ezra^i. In their pioneering work^i^Si^Mi on 
energy flow in the OCS molecule it was shown that LFA 
was capable of highlighting the role of various phase space 
structures including partial barriers and resonance zones. 
The important point is that although motion may be ir- 
regular when viewed over long time scales (several pi- 
coseconds), quantities such as frequencies can be con- 
stant over times corresponding to many vibrational peri- 
ods. This gives rise to long time correlations and the no- 
tion of dynamically significant phase space regiona^L^. 
Since then very little work has been done in utilizing this 
technique for studying energy transport in multimodc 
molecules. Interestingly much more efforts have been 
made in celestial mechanics^ with fewer applications 
in atomic physics^ and molecular systems^. Recently 
Arevalo and Wiggins^^ 9 - have proposed a technique for 
performing time-frequency analysis (LFA) which include 
strongly chaotic trajectories as well by using wavelet 
analysis. Applications to various systems demonstrated 
the advantages of a wavelet based LFA in understanding 
phase space transport. 

A second technique^ - for analyzing the quantum 
states, based on parametric variations, was recently pro- 
posed by us and applied to a 2-mode model Hamiltonian. 
It was shown that variations of the eigenvalues with ap- 
propriate parameters i.e., level velocities were strongly 
correlated to the important classical phase space struc- 
tures organizing the quantum states. In the case of mixed 
and near-integrable systems the magnitudes and relative 
signs of the level velocties clearly indicated^ the dynam- 
ical region of phase space associated with a given eigen- 



state thereby providing a dynamical assignment of the 
states. The technique based on parametric variations is 
dimensionality independent, like the local frequency anal- 
ysis mentioned before, and more importantly basis in- 
variance is guaranteed. Moreover parametric variations 
can be used to identify localization, resulting in devi- 
ations from random matrix theory (RMT) predictions, 
even in strongly chaotic systems^. The full utility of 
parametric techniques is yet to be explored, particularly 
in the context of IVR™, and there do exist detailed semi- 
classical theories in the (near)integrable^ 3 " and chaotic 
situations^. 

Our system of choice DCO is effectively two dimen- 
sional due to the existence of the polyad P. Thus the 
intramolecular dynamics can be investigated in great de- 
tail using the existing methods. However DCO provides 
a stringent test for the dimensionality independent tech- 
niques mentioned above since there is significant chaos in 
the classical phase space for the lowest energies in P = 3 
and 4. This paper is organized as follows: In section 
II the main features of the effective Hamiltonian and its 
classical limit are presented. This is followed, in section 
III, by a dynamical assignment of states in P — 3, 4 us- 
ing the level velocities. It will be shown that a dynamical 
assignment based on level velocities, in conjunction with 
the inverse participation ratios^ (IPR), can be given de- 
spite the presence of significant chaos. In section IV the 
time dependent dynamics of DCO are analyzed using lo- 
cal frequency analysis. This reveals the global picture of 
phase space transport and highlights chaotic trajectories 
undergoing rapid frequency jumps mainly between the 
stretch-stretch 1 : 1 resonance and the DC stretch-bend 
2 : 1 Fermi resonance. Moreover it is seen that chaotic 
trajectories can be trapped in resonances for many vi- 
brational periods. This is related^ to the 'stickiness' of 
the chaotic trajectories and is of importance to the IVR 
dynamics in the molecule. The local frequency analysis 
also highlights the nature of the CO stretch-bend 2 : 1 
Fermi resonance which is seen to play a role in the IVR 
for short times. Section V concludes. 



II. EFFECTIVE HAMILTONIAN FOR DCO 



The spectroscopic Hamiltonian, proposed recently by 
Trollsch and Temps^, was obtained through a fit to the 
experimental DF and SEP spectra. The relevant matrix 
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with j = m,n,b denoting the zeroth order DC stretch, 
the CO stretch, and the bend modes respectively. The 
coupling terms Ai and A2 represent a 1 : 1 and a Darling- 
Dennison 2 : 2 resonance between the stretch modes re- 
spectively. The Fermi resonant couplings k m / n )bb corre- 
spond to a 1 : 2 resonance between the DC (CO) stretch 
and the bend modes. The effective Hamiltonian H com- 
mutes with the polyad Pq = m + n + b/2 and hence block 
diagonal in Pq . A given Pq block has (Pq + 1) (Pq + 2) /2 
states for integer polyads and the states are labelled as 
Pqj in order of decreasing energy. Note that the above 
Hamiltonian is able to approximately reproduce the ex- 
perimental DF spectra for polyad values upto Pq = 5 
since Pq = 6 was the highest polyad considered in the 
fitting procedure. Consequently the current work focuses 
on the polyads Pq = 3,4 as examples. The techniques 
used for the analysis in this work are, however, applicable 
to higher polyads as well. In what follows we will adopt 
a shorthand notation for the two main resonances, 1 : 1 
DC stretch-CO stretch and 1 : 2 DC stretch-DCO bend, 
denoting them as mn and mbb respectively. 

The classical limit Hamiltonian can be obtained using 
the Heisenberg correspondence 



a 7 - <-*■/„•' e 
J j 
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anharmonic resonances. The quantum zero point energy 
Eq has been subtracted from the classical Hamiltonian 
for consistency. The classical system is effectively a two 
degree of freedom system due to the existence of the con- 
served quantity P c = 2(I m + I n ) + I b which is the analog 
of the quantum polyad. As a result the classical dynam- 
ics can be studied in a four dimensional reduced phase 
space. The reduction can be obtained by performing a 
canonical transformation via the generating function: 



between creation-annhilation operators and the action- 
angle variables. The resulting classical Hamiltonian is 

H = ^Ujlj + ^x jk Ijh (3) 

3 jk 

+ Ai(/ m /„) 1/2 cos(0 m - <p n ) + \ 2 I m In cos(2^ TO - 2(j) n ) 
+ -^[kmbblU 2 cos(0 m - 2<j> b ) + k nbb I^ /2 cos(0„ - 2<j) b )] 

with Ai = Ai[l + X[I m + \"I n + A'/'Jfe]. As expected the 
zeroth order Hamiltonian Ho is a function of the actions 
alone and the various couplings manifest themselves as 



2<j) b )J n 



2<j) b )J n + ^ b P c (4) 



The new action-angle variables are related to the old 
ones, J m = I m ,J n = I n ,Pc = 2(J TO + /„) + I b , and 
Tpm = 4>m - 2<f>b,ipn = 4>n ~ 2(f) b ,9 = <j) b . The result- 
ing Hamiltonian2i is independent of the angle variable 9 
which is conjugate to P c and hence an effectively two- 
dimensional system. 

It is possible to perform a Chiriko\*^ analysis of the 
classical Hamiltonian and obtain a resonance template in 
the state spaced. However it is seen that in this case the 
Chirikov analysis, which is based on the zeroth order Hq, 
is not very suitable. This is mainly due to the reduced 
bend anharmonicity of DCO. Despite this limitation use- 
ful qualitative information can be obtained regarding the 
approximate location of the resonance zones and over- 
laps. For instance, such an analysis predicts that the CO 
stretch-bend Fermi resonance (k nbb ) does not show up in 
the state space for most of the polyads. Thus one an- 
ticipates the IVR dynamics in DCO to be dominated by 
the stretch-stretch 1 : 1 and DC stretch-bend 2 : 1 Fermi 
resonances. This observation is in agreement with earlier 
works22*2i. In what follows we will show the resonance 
template as a guideline and the local frequency analy- 
sis of the dynamics establishes the significant deviations 
from the zeroth order picture due to strong couplings. 



III. DYNAMICAL ASSIGNMENTS OF STATES 
IN P = 3 AND P = 4 

Considerable work has been done in the last few 
years to understand and possibly assign highly excited 
states by invoking the classical-quantum correspondence 
principle 26 . Every single method, without exception, is 
based on identifying a dominant classical phase space 
structure that organizes the eigenstates in a given energy 
range. Indeed as DCO is effectively two dimensional due 
to the existence of the polyad it is possible to use any 
one of the several techniques proposed in the literature. 
However our objective in this work is to use techniques 
that do not explicitly rely on determining the exact peri- 
odic orbits and/or visualizations of the phase space and 
wavefunctions. This is especially important for analyzing 
systems with more than two degrees of freedom. 

One such method was recently proposed by us^S and it 
is based on the method of parametric variations. In this 
approach the "level velocities" x a (r) associated with an 
eigenstate |o;(t)) are determined and analyzed. Thus if 
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the effective, spectroscopic, Hamiltonian has the general 
form H = H + ^ . Tj Vj with Vj and Tj denoting the vari- 
ous resonances and the corresponding coupling strengths, 
then the level velocity of a specific eigenstate associated 
with a coupling Vj is obtained as: 

x a (Tj) = ^l = { a \V J \a) (5) 

OTj 

— ^ ] C an Cgn'Vj;nn' — 2 X a (Tj ■; Pj ) (6) 
n.n' Pj 

In the above equation n is the zeroth-order basis and 
c an = (n\a). The level velocity of a state |a) with re- 
spect to Vj can be written as a sum over the polyads 
Pj, corresponding to the resonant perturbation Vj, of 
the partial level velocities x a (rj\ Pj). In the integrable 
limit i.e., with only Vj present it is easy to see that 
^a{Tj) = Xa(i~j;Pj). In the presence of other pertur- 
bations V,^j the level velocity x a (rj) can have contribu- 
tions from many different Pj indicating extensive mixing. 
Thus the level velocity essentially measures the response 
of the quantum eigenvalue to a specific perturbation. 

It was established^ in the earlier work that x a (Tj) is 
correlated to the phase space nature of the eigenstate. 
For integrable, single resonance cases the various proper- 
ties of the level velocities are well known££. In the case 
of mixed phase space systems it is still possible^ to dy- 
namically assign states based on the level velocities since 
the regular regions strongly influence the eigenstates. In 
case of near integrable systems and, to some extent, for 
mixed systems the correlation of the level velocities with 
the phase space can be understood semiclassically since 
it is possible to write^S 

x a (rj) = Tr[V»H] (7) 

* jAtw J didev J w (i,e)w a (i,e) (8) 

where Vj^ is the Weyl symbol of the operator Vj and 
W a is the Wigner function associated with the eigenstate 
| a). For near integrable systems the Wigner function 
condenses on to the invariant torus^ thus explaining the 
observed correlation. 

However for systems with considerable chaos it is not 
a priori clear wether the correlation continues to hold. 
The P = 3,4 cases for DCO considered herein show ex- 
tensive chaos at low energies and hence provide a critical 
test for the level velocity approach. In the extreme limit 
of an ergodic system the Berry- Voros hypothesis states 
that 58 i 59 i 60 

W a (I,6) = CS[E a -H(l,6)} (9) 

In this limit it is easy to see that x a is essentially the 
parametric derivative of the smooth part of the spectral 
staircase function 60 . This term is independent of the na- 
ture of the dynamics and in particular refers to contribu- 
tions from short time orbits. Nevertheless it is well known 



that even in strongly chaotic systems the periodic orbits 
(unstable) play a crucial roleSi. Thus it is expected that 
in such systems the level velocities are influenced by the 
periodic orbits. A more general expectation is that, inde- 
pendent of the nature of the dynamics, closed orbits will 
influence the level velocities. As the purpose of the cur- 
rent work is not to provide a semiclassical expression for 
computing the level velocities, we will briefly outline the 
reasons for such an expectation. The arguments provided 
are not new and appear in an early work by Eckhardt et 
al. where semiclassical matrix elements were determined 
from periodic orbits 6 ^. To begin with define a quantity 

g Vj (E)=TrGVj (10) 

with G being the energy Green's function and Vj is the 
perturbation of interest. Using the formal definition of 
G it is straightforward to show that 

p Vj (E) = ~lmg Vj (E) =^x a {r j )8{E - E a ) (11) 

a 

Thus pvj (E) has poles at the eigenvalues E a and the 
residues are precisely the level velocities x a (rj). The 
quantity gy 6 (E) can be written in the semiclassical limit 
as: 

g V] (E)= J dId6V J w (I,e)G w (I,e;E) (12) 

with Gw being the Wigner transform of the Green's 
function. The above integral can be evaluated via 
stationary phase approximation and, apart from the 
smooth short time contribution, gives rise to an oscillat- 
ing term. The oscillating part can be evaluated for inte- 
grable (Berry- Tabor)S, near-integrable (Almcida-Ullmo- 
Grinberg-Tomsovic)£i, and chaotic (Gutzwiller)S cases. 
In every case the oscillating part depends on the stabil- 
ities of the closed orbits and involves the classical time 
average of Vj over the period of the closed orbits^ 6 . Since 
the level velocities are extracted as residues of the oscil- 
lating part it is reasonable to expect a strong influence of 
the closed orbits and their stabilities. Note that extract- 
ing the velocities from this approach is rather involved. 
Since it is difficult to obtain bounds on the velocities 
using semiclassical techniques in the mixed phase space 
case and the effect of various bifurcations on the level ve- 
locities is not known as of yet we will explore the utility 
of the technique in dynamical assignments from a quali- 
tative perspective. The basis for the assignments will be 
explicitly discussed for the case P = 3 since the argu- 
ments are fairly general. 

In table |U we show the states in polyad P = 3, the 
level velocities with respect to the run and the mbb res- 
onances and the IPR in three different basis. We start 
by noting that the maximum integrable limit velocities 
expected (appendix of Ref. l50l) are 1.5 and 2.8 for the 
x a ,mn and ia^mbb respectively. Moreover since the cou- 
pling constants Ai, k m bb > it is also anticipated, based 
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FIG. 1: State space representation (left panels), Husimi distributions superposed on the surface of section (middle panel), and 
partial level velocities (right panels) are shown for selected states in P = 3. The resonance zones computed using Chirikov 
approximation are also shown in the state space. The mn zone is shown as solid lines and the mbb zone is shown as the dashed 
lines. Three contours of the Husimis are shown for clarity. The states shown are 3i (top row), 34 (middle row), and 3s (bottom 
row). Note that the Husimis of states 3i and 34 are localized around the stable and unstable regions of the classical phase 
space. The Husimi of state 3s is delocalized over the surface of section and the partial velocity data indicates strong mixing. 
The maximum contour value for 3s is one third of the maximum value for state 3i . 



on the earlier work^S, that large positive velocities cor- 
respond to stable (elliptic) regions of the phase space 
and large negative velocities correspond to unstable (hy- 
perbolic) regions of the phase space. Consequently the 
highest energy state 3i is clearly a mn state with the cor- 
responding Husimi localized around a stable region of the 
phase space. The high IPR in the mn basis and the small 
level velocity with respect to the mbb resonance clearly 
support such an assignment. Thus we expect state 3i to 
be localized in the state space as well and associate an 
approximate integrable limit polyad P mn with the state. 
Our expectations are confirmed and Fig^Jshows the state 
space representation of 3i and the corresponding Husimi 
superposed on the surface of section. The approximate 
polyad P mn can be extracted by looking at the partial 
level velocities, shown in Fig^ and one obtains P mn = 3. 
The state 32 has comparable positive velocities both with 



respect to mn and mbb i.e., x m n{^2) ~ imbb{^2)- This 
suggests that the state is influenced by both the reso- 
nances and confirmed by inspecting the state space rep- 
resentation of the state as well as the partial velocities. 
However since the IPR L mn {2t2) > £ m bb(32) and the in- 
tegrable limit mn velocity is expected to be small for this 
state, 32 is assigned as a mn state with P mn = 3. The 
next state 33 can be clearly assigned based on the (par- 
tial) level velocities as well as the IPRs as a mn state 
localized about a stable periodic orbit with P mn = 2. 
State 34 has a negative mn velocity with large IPR in 
the mn basis and hence it is expected to be located in 
the unstable region of the phase space. The Husimi of the 
state shows a typical seperatrix-like structure (cf. Fig^) 
and hence agrees with the velocity analysis. Based on 
the partial velocity data it is straightforward to assign 
an approximate P mn = 3 for 34. Analogous to the state 
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TABLE I: Dynamical assignments for states in P — 3 



State E (cm -1 ) 






Lo 




Lmbb 


Assignment 


3io 


4885 


-0.09 


-1.81 


0.55 


0.56 


0.91 


Pjnbb ; U 


3 9 


4935 


-0.55 


-0.40 


0.65 


0.80 


0.78 


dl 


3 8 


5032 


-0.37 


-0.09 


0.24 


0.33 


0.70 


mbb 


3 7 


5138 


-0.11 


0.45 


0.51 


0.29 


0.46 


*, Pmbb 


3 6 


5177 


-0.02 


1.09 


0.34 


0.33 


0.41 


*) Pmbb )S 


3s 


5289 


-0.31 


0.09 


0.33 


0.60 


0.39 


Pmn 1 


3 4 


5354 


-0.62 


0.23 


0.25 


0.69 


0.22 


Pmn i U 


3 3 


5415 


0.60 


0.13 


0.41 


0.71 


0.34 


Pmn 1) S 


3 2 


5469 


0.22 


0.23 


0.46 


0.81 


0.49 


P m „,nom 


3i 


5644 


1.26 


0.07 


0.54 


0.99 


0.53 


P c s 



a s,u and dl stand for stable, unstable and distorted local respec- 
tively. The approximate polyads P m n = m + n and P mbb = 2m + f>. 
States involved in avoided crossings are denoted by *. 

'The value of P mbb = 6 is deduced from the partial velocity data. 

The value of P mn = 3 is deduced from the partial velocity data. 
See Fig. [T] 



32 it is possible to nominally assign the state 3s as a 
mn state with P mn = 2. State 3io has a large negative 
x m hh and can be clearly associated with the unstable re- 
gion in phase space with P m hb = 6. Again, inspection of 
the Husimi and the partial velocities supports the assign- 
ment. 

We now discuss states that are either difficult to assign 
or unassignable. The particularly strong mixings in these 
states can arise due to avoided crossings and/or strong 
overlap of the various resonances zones. The notoriety of 
multistate avoided crossings in the process of assignments 
are well known and one expects them in sufficiently non- 
integrable systems. State 36 has a large positive mbb ve- 
locity and is expected to be an mbb state localized about 
a stable periodic orbit. However the velocity is about half 
of the maximum expected in the integrable limit and the 
IPRs are about the same in all three basis. The partial 
velocities clearly show that 36 can be assigned as a mbb 
state with P m bb — 6. The state 37 is an interesting state 
since the IPR in the zeroth order basis is larger than the 
IPRs in either the mn or the mbb basis. Indeed the state 
space picture and the Husimis of the state 37 rule out 
a clear mn or mbb assignment. A closer inspection re- 
veals that the states 36 and 37 are involved in an avoided 
crossing with respect to both Ai and k m bb- This explains 
the observed lower than maximum velocity of the state 
36- This quantum mixing, seen in earlier studies2i^&, is 
clearly 2-state and it is possible to demix the states yield- 
ing states \a(f3)) = |3e) ± 1 3y) and analyze the velocities. 
On demixing it is found that x m bb{ot) = 2.2 which is 
closer to the maximum value expected for an mbb state 
in polyad P — 3. The other velocities indicate dominance 
of the mbb coupling and x m bb(P) = —0.65 thus leading 
to an assignment of the state 37 as mbb with P m bb = 4. 

State 38 is a typical example of a difficult state to as- 
sign. In FigQ] the state space representation, Husimi, 



TABLE II: Dynamical assignments for states in P = 4 

State E (cm -1 ) x mn x mbb Lp L mn L mbb Assignment" 

4i 5 6167 -0.67 -0.31 0.80 0.95 0.84 dl 

4i4 6392 -0.30 -3.38 0.22 0.25 0.76 P m bb,u 

4i 3 6448 -0.30 -0.98 0.36 0.40 0.76 P mbb 

4 12 6593 -0.30 0.83 0.22 0.27 0.81 P mbb 

4 n 6715 0.07 -1.35 0.54 0.50 0.68 P mbb -2,u 

4i 6773 -0.28 3.35 0.32 0.32 0.72 P mb b,s 

4 9 6809 -0.59 -0.47 0.38 0.70 0.47 dl 

4 8 6875 -0.37 0.13 0.13 0.29 0.34 mixed 

4 7 6965 0.15 0.09 0.47 0.30 0.35 *,P mn -2 

4 6 7017 -0.08 1.18 0.30 0.27 0.27 *,P mbb -2,s 

4 5 7101 -0.01 0.32 0.36 0.57 0.40 mixed 

4 4 7168 -0.81 0.19 0.27 0.78 0.26 P mn ,u 

4 3 7241 1.03 -0.06 0.33 0.61 0.41 P mn -l,s 

4 2 7287 0.64 0.41 0.23 0.65 0.29 P mn 

4i 7470 1.81 0.07 0.49 0.99 0.47 P mn ,s 

"The notations are same as in table |I] The values P mn = 
4, P m bb = 8 are deduced from the partial velocity data. See Fig. 121 

and the partial velocity data is shown for 3s- The veloc- 
ity data indicates this to be a mn state whereas the IPR 
data supports a mbb assignment. Note that this state 
has the lowest IPR in the zeroth order basis amongst the 
set of states in P = 3. The particularly high value of the 
IPR in the mbb basis suggests that it is quite possible 
for this state to have a small velocity even in the single, 
mbb only integrable case. Indeed on computing the in- 
tegrable limit one obtains imbbi^s) ~ —0.13. Thus it is 
possible to assign the state 3s as a mbb state. Note that 
the partial velocity data suggests P mn = 4 but the state 
space representation exhibits P m bb = 6. The difficulty, 
in part, arises from the fact that it is quite possible for a 
state to have a velocity close to zero and yet be strongly 
influenced by that resonance^. A similar situation arises 
for the state 3g in that the level velocities are very simi- 
lar. In this case, however, the IPRs are very similar and 
high in all three basis. This is a typical signature of a 
distorted local mode state. In this instance the state is 
influenced slightly by both the resonances leading to the 
similar velocities and IPRs. 

The states for P = 4 can be assigned in a similar fash- 
ion and the assignments are shown in table [H] In par- 
ticular, all the states except 4§ and 4s can be assigned 
in analogy with the states in P = 3. Thus the states 4s 
and 45 are classified as mixed. In Fig[21 states with large 
positive (4io) and large negative (^14) mbb velocities are 
shown. Note the localization of the Husimis in the stable 
and unstable regions of the phase space. The state space, 
Husimi and partial velocity data for the mixed state 4g 
are shown in Fig|2]as an example. Extensive delocaliza- 
tion is seen for this state in both the state space and the 
phase space. The state 4g also happens to have the low- 
est IPR among the states in P = 4. We note that for 



6167 


-0.67 


-0.31 


0.80 


0.95 


0.84 


dl 


6392 


-0.30 


-3.38 


0.22 


0.25 


0.76 


Pmbb 5 ^ 


6448 


-0.30 


-0.98 


0.36 


0.40 


0.76 




6593 


-0.30 


0.83 


0.22 


0.27 


0.81 


P"mbb 


6715 


0.07 


-1.35 


0.54 


0.50 


0.68 


■Pmbb 2^ 


6773 


-0.28 


3.35 


0.32 


0.32 


0.72 


Pmbb 3 S 


6809 


-0.59 


-0.47 


0.38 


0.70 


0.47 


dl 


6875 


-0.37 


0.13 


0.13 


0.29 


0.34 


mixed 


6965 


0.15 


0.09 


0.47 


0.30 


0.35 


* 5 Pmn 2 


7017 


-0.08 


1.18 


0.30 


0.27 


0.27 


*j Pmbb 2, ; 


7101 


-0.01 


0.32 


0.36 


0.57 


0.40 


mixed 


7168 


-0.81 


0.19 


0.27 


0.78 


0.26 


Pmn t U 


7241 


1.03 


-0.06 


0.33 


0.61 


0.41 


Pmn lj S 


7287 


0.64 


0.41 


0.23 


0.65 


0.29 


P 

- 1 mn 


7470 


1.81 


0.07 


0.49 


0.99 


0.47 


Pmn •> S 



7 




FIG. 2: Same as in Fig. Q for three states in P — 4. The states shown are 4io (top row), 4i4 (middle row), and 4g (bottom 
row). Note that the Husimis of states 4io and 4 14 are localized around the stable and unstable regions of the classical phase 
space. The Husimi of state 4s is delocalized over the surface of section and the partial velocity data indicates strong mixing. 



P = 4 the maximum expected velocities for the mn and 
the mbb are 2.0 and 4.4 respectively. 

A general observation regarding the assignments is 
that the lower energy states of the polyads are domi- 
nated by mbb whereas the higher end states are influ- 
enced more by the mn resonance. States in the middle 
of the energy regions exhibit transitional character. Simi- 
lar observations were made in the dynamical assignments 
of the states in P = 8 by Jung et at—. Note that the 
assignments made in this work are essentially based on 
the level velocities and the IPRs. Thus our assignment 
does not invoke any explicit quantum numbers based on 
periodic orbits. However it is clear that attaching an 
approximate polyad and noting the signs on the level ve- 
locities is providing dynamical information. For instance 
state 4i having a large i mn > implies that the state 
is influenced by a stable periodic orbit of the ran type. 
In particular this state would be considered as a C class 
state with the assignment (I = 4, t = 0) in the notation 
of the previous work- 1 . The approximate polyad P mn 
can be identified with the longitudinal quantum number 



I. A crucial point regarding the assignments using level 
velocities is worth reiterating at this juncture. The ad- 
vantage of using the level velocities is based on strong 
correlation of the velocities with the phase space nature 
of the eigenstates. Such dynamical informations cannot 
be obtained by using the IPRs alone or the various rep- 
resentations of the eigenstates. For example, in P = 3 
the IPRs imn (33) ~ (34) and one would correctly 
infer that both states are of mn type. However with 
Xmniii) ~ —x mn (^3) it is possible to provide a finer dis- 
tinction between the two states ; 33 is associated with a 
stable region of the phase space whereas 34 is located in 
the unstable region of the phase space. 



IV. LOCAL FREQUENCY ANALYSIS: PHASE 
SPACE DIFFUSION AND TRANSPORT 

As mentioned before the Chirikov analysis of DCO res- 
onances based on the zeroth order classical Hq(I) is not 
very accurate due to strong resonant couplings of the 



8 



D.6 




t(ps) 



FIG. 3: Survival probabilities for the ZOBS (a) |030) and (b) 
1 040). Solid line is the exact Hamiltonian and dashed lines 
are with k n tb — 0. 



various modes. The surface of sections indicate that for 
polyads P = 3,4 there is considerable stochasticity even 
at the lowest energy. This implies that the various reso- 
nances overlap strongly for these polyads and most tra- 
jectories are chaotic. Despite the significant amount of 
chaos in the phase space the Husimis of many eigenstates 
were seen to be rather localized. Such observations have 
been made earlier in other studies and it is well known 
that a largely chaotic phase space can still contain struc- 
tures that act as partial barriers to phase space transport. 
The partial barriers can arise due to cantori as well as 
broken seperatrices and existence of such structures can 
lead to eigenstate localization. Moreover even in the case 
of strongly chaotic systems eigenstates can get scarred by 
periodic orbits leading to localization. In recent years it 
has been understood that such partial barriers and peri- 
odic orbit scarring lead to statistically significant devia- 
tions from the predictions of RMT. Infact the level veloc- 
ities, studied in the previous section, already encode such 
localization information. From the IVR perspective, sev- 
eral studies^ have now established that the energy flow 
process is hierarchical and that a hierarchical local ran- 
dom matrix (HLRM) approach is more appropriate to 
understand IVR. The important factor is the local den- 
sity of coupled states and not the total density of states 
at a given energy which leads, amongst other observa- 
tions, to a power law decay of the survival probability on 
intermediate timescales^i. 

Based on the detailed work on classical phase space 
transport^ it is tempting to conjecture that the hier- 
archical nature of IVR is intimately connected to the 
various partial barriers in the phase space. The im- 



portant wor 



41.42 



by Davis and Martens et al on OCS 



seems to support the conjecture. However, except for 



earlier works on OCS, there have not been many efforts 
to identify partial barriers for transport in the context 
of IVR in molecules. The importance of such an un- 
dertaking lies in the fact that it might prove crucial in 
controlling IVR as well. In DCO it has been established 
that the unimolecular decay is IVR limitedl£i2il. More- 
over, previous studiesi&22*21 on DCO suggest that the 
Fermi resonance between the CO stretch and DCO bend 
is unimportant. The approximate Chirikov analysis also 
indicates the muted role of the nbb Fermi resonance. In 
FigGDthe survival probabilities for the zeroth order bright 
states 1 030) and 1 040) in the absence of the nbb Fermi res- 
onance are shown and compared to the full system. It is 
clear from the figure that although the initial decay from 
the bright states is controlled by the 1 : 1 (for about 100 
fs) the nbb Fermi resonance starts to play a role around 
200 fs. Is it the case that the nbb resonance does play 
a role, albeit in a secondary fashion, in the short time 
tier structure? Why is the initial decay from 1 040) ap- 
parently complete by about 150 fs as opposed to those 
of 1 030) and |050) ? In other words what makes the 1 : 1 
so effective for the bright state |040) ? The answers to 
these questions are not straightforward and in this work 
we will try to provide some approximate answers. We 
believe that the LFA technique is capable of addressing 
some of these crucial issues. Ahead a brief introduction 
is given to the technique of LFA and we refer the reader 
to the work by Martens et a& and Arevalo et a£L2£ for 
more details. 

The classical limit of the spectroscopic Hamiltoni- 
ans arc nonlinear, multiresonant Hamiltonians which are 
functions of the zeroth order actions I and the angles 
(p. The classical limit Hamiltonians have a typical form 
H(I,4>) = H Q (I) + ^jVj(l,(j>) where the various reso- 
nant perturbations are denoted by Vj . In the absence of 
the perturbations the system is described by Hq which 
is classically integrable. In such cases the various nonlin- 
ear frequencies = dHo/dl are time-independent and 
one can define a frequency map I — * fi(I) which, if invert- 
ible, can be used to characterize the phase space. The fre- 
quency map can be used to characterize the phase space, 
thanks to the KAM theorem, for near-integrable systems 
as well. For a system like DCO there are three frequen- 
cies f2 = (O m , fin, rib) and it is useful to visualize the 
dynamics in the frequency ratio space. However for non- 
quasiperiodic trajectories, associated with the full nonin- 
tegrable Hamiltonian, the various frequencies 51(1) vary 
with time since increasing perturbation strengths lead to 
the breakdown of KAM tori and widespread chaos. De- 
spite the time dependence of the frequencies it has been 
observed in many systems that the frequencies can re- 
main approximately constant over times corresponding to 
many vibrational periods. Thus even non-quasiperiodic 
trajectories can spend an appreciable amount of time in 
one part of the phase space before passing on to another 
region of the phase space. This in turn implies existence 
of long time correlations and gives rise to the concept 
of dynamically significant regions of phase space. The 
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concept of local frequencies i.e., time- varying frequen- 
cies allows one to identify such dynamically significant 
regions of phase space. 

Two main approaches are available to extract the lo- 
cal frequencies from a classical trajectory. The first 
one proposed by Martens et aiai, and a similar one by 
Laskar^, uses Fourier transforms of short time segments 
of long trajectories. We are interested in the more recent 
work4£f^ by Arevalo and Wiggins wherein the local fre- 
quencies are obtained by a continuous wavelet transform 
of the trajectories. In this approach a general time de- 
pendent function f(t) is expressed in terms of the basis 
functions constructed as translations and dilations of a 
mother wavelet^ 



1pa,b(t) 



-1/2 



lb 



t 



(13) 



The coefficients of this expansion are given by the wavelet 
transform of f(t), defined as 



L^f(a, b) = a 



-1/2 



dt 



(14) 



for a > and b real. The wavelet transform gives the lo- 
cal frequency of f(t) over a small interval of time around 
t = b and inverse of the scale factor a is proportional to 
the frequency. A distinct advantage of the wavelet ap- 
proach is that the time window automatically narrows for 
high frequency and widens for low frequency. The value 
of the local frequency at a given time t = b is extracted^! 
by determining the scale a which maximizes the modu- 
lus of the wavelet transform. Throughout this work, as in 
Arevalo and Wiggins, we will use the Morlet-Grossman 
mother wavelet 



m = 



i 



erV27r 



e 2iriAt g-r /2a 2 



(15) 



with A = 1 and a — 2. The parameters A, a can be tuned 
to improve the resolution. Within the wavelet approach 
the local frequencies can be obtained even for chaotic 
trajectories. For a strongly chaotic trajectory the fre- 
quencies will vary considerably and one can measure the 
"diffusion" of the k th frequency by evaluating 



d k (T) = ^ I ,//|<>,(/) -<>| 



(16) 



where Cl is the mean frequency and T is the total time to 
which the trajectory is propagated. This measure of dif- 
fusion characterizes the unstable regions of phase space 
and thus low diffusion implies either quasiperiodic motion 
or extensive trapping near resonance zones (see below). 
For our application to DCO we choose initial conditions 
for the trajectories in the following fashion. A uniform 
grid is chosen in the state (zeroth order action) space 
{Im,In) with lb being fixed by the polyad conservation. 
A particular choice of the angles {4> m ,4>rn4'b) then pro- 
vides initial conditions corresponding to a specific slice 
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FIG. 4: The three frequency ratios for DCO for various tra- 
jectories at energy E = 5145 cm -1 . Total time of propagation 
T = 10 ps. Qm/iln is shown in blue, £l m /£lb is shown in red 
and Q.n/Q.b is in green. See text for details on each trajectory. 



of the phase space. The numerically integrated trajec- 
tories over certain time interval [0, T] are expressed as 
z k (t) = ^2/ fc (t)e^ fc(t) for k = m,n,b. The three local 
frequencies f2fc(i) are extracted by computing the max- 
ima of the modulus of the wavelet transforms of Zk(t) i.e.. 
max a |L^,Zfc(a, b)\. Note that the maxima are obtained 
with a fractional precision of 10~ 8 . Details regarding the 
numerical procedure to extract the maxima can be found 
in the work of Arevalo et atH. 

A large number of trajectories for both P = 3 and 
P — 4 at various energies were analysed and many dis- 
tinct dynamical behaviours were repeatedly seen. In 
Fig0|we show, as examples of such dynamical patterns, 
six trajectories whose frequency ratios were computed to 
about T — 10 ps. Note that 10 ps is a fairly long time 
and corresponds to about 500 CO vibrational periods. 
All of the six trajectories are at an energy E = 5145 
cm -1 which is energetically close to the middle of the 
polyad P = 3. The phase space structure is shown in 
Fig|S]for reference and one can clearly see large regions 
of stochasticity. Figure (a) shows a typical quasiperiodic, 
nonresonant trajectory which in this case represents a 
distorted local mode. The corresponding trajectory is 
shown in the surface of section (circles in Fig0 as well. 
In Fig^3 a 1 : 1 resonant regular trajectory correspond- 
ing to the small island on the surface of section is shown. 
Notice the locking of the relevant frequency ratio i.e., 
Qm/^n ~ 1 for the entire duration. In FigQJi a trajec- 
tory with initial conditions very close to a periodic orbit 
of the full system is shown. All the frequency ratios are 
locked perfectly for about 7 ps and later show large diffu- 
sion. Interestingly in this case the nbb locking fl n /Qb = 2 
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FIG. 5: Classical Poincare surface of section at E = 5145 
cm -1 . The variables are J m = I m and tp m — <p m — 2(f>b- 
Selected trajectories (cf. Fig0J are highlighted with filled 
circles (a), triangles (d), and squares (e). The highlighted 
trajectories are plotted for about 5 ps. 



is clearly seen and thus any other initial condition close 
to this would certainly involve the nbb resonance to some 
extent. In Fig^Ji we show an example of a trajectory ex- 
hibiting resonance trapping. This particular trajectory 
is trapped in the 1 : 1 resonance for nearly 5 ps and then 
drifts away. The corresponding trajectory is plotted in 
the surface of section (triangles in FigEJ and it is seen 
that this trapping is due to the stickiness of the small 1 : 1 
island. The trajectory is chaotic for longer times but the 
observed stickiness suggests that even tiny regions of reg- 
ularity in the phase space can sufficiently influence the 
dynamics. Another example of a trapped trajectory is 
shown in Fig0Jj with the corresponding behaviour on the 
surface of section shown in Fig|5] (squares). In this case 
the trajectory starts out with frequency locked in the 
1 : 1 but drifts away quickly and gets trapped into the 
Qm/^b ~ 2 ratio for almost 6 ps. Finally in Fig0f we 
show a typical chaotic trajectory with large diffusion and 
no apparent long time trapping in any of the fundamental 
resonances. Two points are worth reiterating in the con- 
text of the examples shown in FigQJand Fig[5] The first 
is that all these various, and other, dynamical behaviors 
are found at the same energy indicating the rich struc- 
ture inherent to the classical phase space. This energy 
is by no means special and indeed similar observations 
hold for other energies in the polyad P = 3 as well as dif- 
ferent polyads. Secondly, the overall phase space might 
look chaotic on long time scales but the regular regions, 
even those with phase space area < h, lead to significant 
stickinesa^SiS and hence long time trappings in one or 
the other significant resonance zone. This observation, 
we emphasize again, is not new and is well known in the 




t(ps) 

FIG. 6: A trajectory for P = 4 exhibiting rapid frequency 
jumps. For clarity the ratios S7 m /f2(, (solid line with points) 
and Q.n/£lb (dashed) are shown. Note that the jumps (in- 
dicated by arrows) take the trajectory from being near mbb 
resonance to being near the mn resonance. Only a 2.5 ps por- 
tion of the 15 ps trajectory is shown. The jumps occur over 
the entire time span. 



nonlinear dynamics community. However in this work we 
are demonstrating this effect for a realistic spectroscopic 
Hamiltonian with important consequences for IVR in the 
molecule. 

A further dynamical effect that occurs in DCO has to 
do with rapid jumps in the frequency ratio space. An 
example for P = 4 is shown in the Fig|fj] where the ra- 
tio f2 m /f2{, is undergoing rapid jumps between otherwise 
smooth variations. Note that the ratio is fairly 

smooth and does not exhibit the rapid jumps. The jumps 
in the frequency ratio happen from the trajectory being 
near the mbb resonance to the trajectory being near the 
mn resonance. Such rapid jumps have been observed by 
Arevalo in OCS as well as the planar 3-body problem 49 . 
In this instance the jumps seem to correspond to hete- 
roclinic orbits that connect the mbb and the mn reso- 
nance zones. However confirming this involves comput- 
ing invariant manifolds and their intersections which is 
far from easy. Nevertheless, such rapid resonance tran- 
sitions could have important implications for the nature 
of IVR in DCO. 

To indicate the strongly coupled nature of the DCO 
dynamics we show the total diffusion for P = 4 over the 
state space in FigH In FigE^ the diffusion is shown for 
the phase space slice cf> m — cf> n = <\> b = 0. For refer- 
ence Figs[7|;,d show the diffusion data for integrable mn 
and mbb subsystems. It is easy to see the mn and mbb 
resonance zones in the integrable single resonance cases. 
Comparing to FigE^, it is clear that the resonances have 
overlapped strongly. Further in FigEt> we show the diffu- 
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FIG. 7: Total diffusion (ps^" 1 ) in the state space for P = 4 computed for T = 10 ps. In (a) and (b) the phase space slices are 
(0, 0, 0) and (0, tt, 0) respectively. The integrable subsystems are shown for reference in (c) mn and (d) mbb respectively. In 
(c) and (d) the resonance zones can be seen clearly. In (a) and (d) there are no isolated resonances and the system is strongly 
coupled. 



sion for the phase space slice <p m — = 4>b, 4> n = tt. Note 
that the diffusions for the two slices of phase space are 
dramatically different. This implies that it is not possi- 
ble to understand the intramolecular dynamics of DCO 
based on Ho alone. A key observation here is that the 
diffusions in state space show significant heterogeneity. 
For instance in Fig[7Ji,b we have regions of similar dif- 
fusion strengths seperated by a region of larger diffusion 
strengths. Thus if one imagines going along a specific mn 
or mbb polyad in state space then alternating diffusion 
strengths are encountered. Similarly along a specific en- 
ergy contour the diffusion values can vary considerably. 
This leads to trappings of the classical trajectory and 
hence bottlenecks to energy flow. For instance in Fig[7|3 
the region of low diffusion around (I TO , I n ) » (1.0, 0.5) is 
due to the existence of a periodic orbit of the full system. 
In the semiclassical limit the heterogeneous nature of the 
state space diffusion can have interesting consequences 
for eigenstate localizations. 

Temps et al. discuss^ the differences between choos- 



ing a CO-stretch bright state (experimental) versus the 
DCO bend bright states. It was shown that the DCO 
bend states undergo much faster unimolecular decay. In 
order to gain dynamical insights into their findings the 
various bright states were analyzed from the LFA per- 
spective. Several initial conditions corresponding to the 
bright states i.e., E = Ezobs were generated with the 
actions (I m , I n , If,) corresponding to the bright state of 
interest and randomly choosing the angle (<p m ,4 > n,<Pb) 
variables. The various dynamical quantities of interest 
are then averaged over the angles. This is done for two 
main reasons. First, although the correspondence be- 
tween actions and quantum numbers is well known a sim- 
ilar correspondence for angles is not very clear. In the 
deep semiclassical limit one expects a correspondence but 
the problem at hand is far from such a limit. Second ob- 
servation is that the quantum dynamics associated with 
a particular ZOBS, semiclassically speaking, would nat- 
urally involve a family of classical trajectories which are 
modeled here with the random angle distributions. The 
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FIG. 8: Average crossings per unit time (a), average longest 
locking time (b), total time near a specific resonance (c) and 
diffusion of the frequencies (d) are shown in this figure for 
six ZOBS. A, B, C, D, E, F correspond to the bright states 
|030), |040), |050>, |006), |008>, and |0010> respectively. Aver- 
aging is done over 1000 trajectories and T = 10 ps. In (a), 
(b), and (c) the filled circles and triangles correspond to the 
mn and mbb resonances respectively while the open symbols 
are for the case with k n bb — 0. See text for details. 



fact that the diffusion of the frequency ratios are strongly 
dependent on the choice of the phase space slice further 
support the procedure of angle averaging. The quan- 
tities of interest are number of crossings per unit time 
of the two main resonances (denoted (c mn , c m bb) /T), to- 
tal time spent near the resonances (denoted (r mra , T m u>)), 
longest locking time in a given resonance (denoted (r m )), 
and the diffusions d(T). In this study we have averaged 
over 1000 trajectories with the angles chosen such that 
E = Ezobs ±1.0 cm -1 and each trajectory was propa- 
gated to T = 10 ps. 

In FiglHl we show the various dynamical quantities of 
interest for three of the experimental CO stretch bright 
states and three of the bend only bright states corre- 
sponding to the polyads P = 3, 4, 5. In Fig|S^ the num- 
ber of resonance crossings per unit time are shown. It is 
immediately clear that the mn resonance plays a domi- 
nant role for the CO stretch bright states and that the 
mbb resonance is key for the bend bright states. Notice 
that for the |008) and the |0010) bright states the 'flux' 
across mn is almost zero. On the other hand despite 
the dominance of mn flux for the CO stetch ZOBS, the 
mbb flux is not negligible. This already provides a hint 
to the drastically different IVR dynamics from the two 
classes of ZOBS with |003) being an exception. Fig|S|3 
shows the longest locking time, on average, for the dif- 
ferent bright states. Interestingly even though the mn 
flux is larger than the mbb flux for the CO stretch states 



the r m 's show the opposite trend. The total time spent 
near a particular resonance as shown in FigEb a lso in- 
dicates that the IVR from CO stretch states is strongly 
effected by the mbb resonance. The diffusions in FigQJl 
show marginal differences. A clear picture emerges for 
the IVR dynamics associated with the CO stretch bright 
states. Initially the mn resonance transfers energy to the 
DC stretch (dissociation mode) but prior to a substantial 
buildup of excitation in the DC stretch the energy is "si- 
phoned off" into the bend mode via the mbb resonance. 
This process happens for a fairly long time leading to 
a highly nonstatistical nature of the dynamics and uni- 
molecular decay of DCO. On the other hand for the bend 
bright states corresponding to P = 4, 5 the initial IVR 
dynamics is dominated by the mbb resonance. However 
the mn resonance is unable to compete and this results 
in unhindered transfer of energy to the dissociative mode 
via the mbb coupling. Exception to the nature of bend 
ZOBS is clearly seen in the case of |006) wherein the mn 
resonance does compete and leads to transfer of energy 
from the dissociative mode to the CO stretch mode. An 
important observation is that the diffusion of frequencies 
computed here is the lowest for the ZOBS |0010) which 
undergoes the fastest unimolecular decay. This is not a 
contradiction since large difusion corresponds to many 
hops and consequent locking leading to incomplete flow 
of energy into the dissociative mode. For example the 
largest diffusion is predicted for the ZOBS |006) which 
is consistent with the fact that the times spent near the 
two main resonances are similar as compared to the other 
states. Thus |006) should decay the slowest which is con- 
firmed by computing the associated survival probability 
(cf. figure 3(b) in Ref . I20I) . It is interesting to observe 
that the diffusions for the CO stretch bright states are 
not very different whereas in the case of the bend bright 
states there is a large variation. This demonstrates the 
mode-specificity of the intramolecular dynamics of DCO. 



In Fig|Slthe effect of setting k n bb — is shown. It is 
clear that the largest effect of the nbb Fermi resonance 
is for the CO stretch bright states |030) and |040). In 
the absence of nbb the mn flux is decreased and the mbb 
flux is increased with larger time being spent, on the 
average, near the mbb Fermi resonance. This indicates 
that the CO stretch bright states would lead to smaller 
rates of unimolecular decay in the absence of the nbb 
resonance. This is true for the bend bright states as well 
although the effect is very small. Thus the role of nbb 
Fermi resonance is not negligible for the lower polyads. 
As a final observation we point out that amongst the CO 
stretch bright states there is a sharp increase, on going 
from 1 030) to the |040), in the average residence time and 
the average longest locking time for the mn resonance. 
This perhaps explains the efficiency of the mn resonance 
in P = 4 but at this stage it is conjectural. 
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V. CONCLUSIONS 

A long term goal of our research is to come up with 
techniques that would allow fairly detailed classical- 
quantum correspondence studies of highly excited states 
in systems with three degrees of freedom or more. To this 
end two dimensionality independent techiniques, para- 
metric variations and local frequency analysis, have been 
employed in this study. We have demonstrated the utility 
of both the techniques for understanding the intramolec- 
ular dynamics of DCO. The successes of both the tech- 
niques, in dynamical assignments of states and establish- 
ing the nonstatistical nature of the dynamics, given the 
strongly coupled nature of the system (coupling strengths 
comparable to the mean level spacing) is quite encourag- 
ing. We conclude with brief comments on the techniques 
and future outlook. 

The level velocity approach, used to dynamically as- 
sign the eigenstates, is certainly capable of identifying 
strongly localized states. The close connections between 
the level velocity of an eigenstate and its corresponding 
phase space nature makes this approach attractive. Pre- 
vious application^ to a mixed 2-mode system demon- 
strated the utility of the level velocities for dynamical 
assignments of the highly excited states. In this work 
we have shown that the method is useful even in the 
presence of significant chaos. Although avoided crossings 
lead to difficulties, and this is not surprising, it is debat- 
able wether the resulting mixed states can be or should 
be assigned. Demixing the states invovled in avoided 
crossings and then analyzing them can be done within 
this approach as in any other approach. However the 
procedure of demixing, especially in multistate or broad 
avoided crossings, should not be considered as an assign- 
ment. This work shows that analyzing the level veloc- 
ities in conjunction with measures like IPRs can be a 
fairly useful way of dynamically assigning highly excited 
states. At this stage the effects of various bifurcations on 
the level velocity spectrum is not known. Interestingly in 
the case of DCO, despite substantial change in the phase 
space structure with changing energies, the level velocity 
corresponding to a state localized about stable regions 
is very close to the classical estimate. This has been 
observed before^ in the case of HCP as well. This sug- 
gests a very detailed classical-quantum correspondence 
and needs to be studied further. In order to make this ap- 
proach more useful it is important to obtain some bounds 
on the level velocities in the mixed and strongly chaotic 
situations. Thus a semiclassical theory for the velocities 
in the chaotic regimes needs to be developed and work is 
in progress in our group along these lines. Note that the 
distribution of the velocities has been found to be gaus- 
sian in the strong chaos limitSi. The general technique 
of parametric variations has a wide range of applications 
and level curvatures have also been studied 7 ^. It remains 
to be seen if the level curvatures have any utility in un- 
derstanding highly excited states and IVR. For example 
conductance of disordered systems have been related to 



level curvatures£*Z2i£i and given the tantalizing similar- 
ities between the example mentioned and IVR it seems 
worthwhile to investigate the utility of level curvatures 
in molecular systems. 

The LEA technique is able to provide detailed phase 
space information in a strongly coupled system like DCO. 
In particular the time-frequency analysis revealed the ex- 
istence of trajectories trapped near important resonance 
zones. In an effective spectroscopic Hamiltonian the key 
resonances are easy to expect since the fit already builds 
in the main resonant perturbations but LFA is able to 
provide information on the competing roles of these var- 
ious resonances. An example is the clear differences be- 
tween the CO-stretch bright state and DCO bend bright 
state dynamics for DCO. LFA is also able to highlight 
the role of weak resonances like the nbb Fermi resonance. 
As shown in this work and in the previous works the 
wavelet based approach to LFA is useful in analyzing 
effective Hamiltonians and various dynamical quantities 
can be determined to understand IVR. LFA clearly shows 
that despite significant chaos the DCO phase space is far 
from the RMT regime. In this paper it is shown that even 
for a strongly coupled system like DCO the state space 
diffusion has a lot of structure with regions of high and 
low diffusions interspersed thorughout. The diffusion of 
frequencies provides clear insights into the dynamics at 
a particular polyad and/or a specific energy. This infor- 
mation can then be used to determine the nature of the 
'random walk' in action space for strongly coupled sys- 
tems and thus provide more insights into certain recent 
conjectures^ regarding the intermediate time behavior of 
survival probabilities. Although not shown here most of 
the trajectories exhibit across-resonance diffusion which 
is quite fast and implicates higher order secondary res- 
onances as well 4 ^* 7 ^. The rapid frequency jumps which 
perhaps correspond to heteroclinic orbits joining unsta- 
ble periodic orbits have not been investigated in detail. 
Surely the rapid jumps from one resonance to another 
have important consequences for IVR and needs further 
study. We also note that there is not much of an Arnol'd 
web for DCO for P = 3,4 in contrast 47 to H 2 0. DCO 
being effectively two dimensional implies that there is 
no transport along resonance lines and thus many of 
the interesting conjectures put forward^ 1 ^ by Martens, 
Davis, and Ezra remain to be tested. Currently work is 
in progress for true three degree of freedom systems. 

As a final note we mention that both techniques used 
in the current study probe localiztion characteristics. 
The parametric response is probing eigenstate localiza- 
tion and the LFA is probing localization of phase space 
transport due to partial barriers and broken seperatri- 
ces. Localization explicitly implies nonstatistical behav- 
ior and thus deviations from RMT and incomplete IVR in 
the system. Evidence existsZiZi that the quantized ver- 
sions of the partial barriers can be more restrictive than 
the classical analogs. A interesting link between these 
two approaches is provided by the recently proposed 
intensity-velocity correlator^iiSiS. Detailed studies of 
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true three degrees of freedom system via the intensity- 
velocity correlator can lead to a better understanding of 
the IVR dynamics. 
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